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Graphene-based materials show promise for spintronic applications due to their potentially large 
spin coherence length. On the other hand, because of their small intrinsic spin-orbit interaction, 
an external magnetic source is desirable in order to perform spin manipulation. Because of the 
flat nature of graphene, the proximity interaction with a ferromagnetic insulator (FI) surface seems 
a natural way to introduce magnetic properties into graphene. Exploiting the peculiar electronic 
properties of bilayer graphene coupled with FIs, we show that it is possible to devise very efficient 
gate-tunable spin-rotators and spin-filters in a parameter regime of experimental feasibility. We also 
analyze the composition of the two spintronic building blocks in a spin-field-effect transistor. 

PACS numbers: 72.25.-b, 72.80.Vp, 85.75.-d, 85.75.Mm, 85.70.-w 



Graphene with its high mobility-^ and potentially long 
spin lifetimes, is an attractive material for spintronics. 
In particular, spin relaxation lengths on the order of mi- 
crometers have been observed^, together with spin re- 
laxation times of hundreds of picoseconds, which are 
still believed to be limited by extrinsic impurities^. 
More recent experiments reported the measurement of 
a spin lifetime up to 1 ns in graphene and even of sev- 
eral nanoseconds in bilayer graphene (BG^. Moreover, 
tunnel-injection of spin into graphene has been recently 
achieved using Co ferromagnets, with the observation of 
the largest non-local magnetoresistance of any material^. 
Graphene quantum dots have been also identified as an 
ideal host for spin qubits^. 

The reason for such favorable spin properties is the 
small spin-orbit coupling (SOC) and the weak hypcrfinc 
interaction with the underlying nuclear spin syste m 10 ! 11 . 
The SOC in single layer graphene has been predicted to 
be on the order of 1CP 3 -10 -2 meV— ~— . On the other 
hand, this weak SOC constitutes a severe limitation for 
spin manipulation in conventional spintronics devices like 
the Datta-Das spin-field-effect transistor (SFET) 15 . 

An alternative strategy is offered by contacting 
graphene with a ferromagnetic insulator (FI) , giving rise 
to an exchange proximity interaction (EPI)^&~— . EPI re- 
sults from the Coulomb exchange interaction between 7r 
states in graphene and the magnetic ions on the FI sur- 
face. Ideally, the EPI — being short-ranged — affects only 
a graphene layer in direct contact with the FI and acts 
like an effective Zeeman field, superimposed on the orig- 
inal BG Hamiltoniar>i£. 

Here, we theoretically study transport through BG in 
a double gate configuration, on which a FI is used as 
spacer between the upper (U) layer and the top gate, 
giving rise to EPI, as shown in Fig. Ufa). The gate bias 
A is used to impose a semiconducting gap and to localize 
the low-energy region of the conduction band on the U or 
the lower (L) layer (depending on the sign of A). Indeed, 
a tunable semiconducting gap up to 250 meV with the 
application of a gate bias has been demonstrated^— . 
Consequently, it is possible to electrically control the ef- 



fective Zeeman field for electrons in the conduction band, 
turning ON or OFF the device. We show, in particu- 
lar, that the device can act either as a spin-rotator (SR) 
or as a spin-filter (SF). Finally, we propose and analyze 
the combination of these two spintronic building-blocks 
within a SFET. 



I. SETUP AND SCATTERING PROBLEM 



The Dirac Hamiltonian describing BG near the K point 
H = -^-t x + v f {a x p x + (T y T z Py) + -^-(a z + (Tq)t x , (1) 



with a and r the Pauli matrices for the sublattice 
(A,B) and layer (U,L) degrees of freedom, A the po- 
tential energy difference between the U and L plane, 
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FIG. 1: (Color online) (a) Setup of the double gate BG with 
FI used as a spacer between the U layer and the top gate, (b) 
Dispersion curve of biased BG, showing the typical Mexican 
hat behavior, subject to EPI. A gate-bias A = —150 meV is 
imposed between the BG layers. The full blue lines are /*•- 
spin polarized bands (concord to the FI magnetization axis), 
while the dashed red lines are ^-spin polarized. The dotted 
lines represent the normal non-magnetic BG dispersion, as 
used for the LS and RS leads. 



t± = 0.39 eV the inter-layer hopping parameter—, and 
Vf ~ 10 6 m/s 24 . The Hamiltonian acts on the the spinor 
* = (X-B'i XA', XA: Xb), with A and -B on the U layer, A' 
and B' on the L one. 

When the U layer is placed in direct contact with the 
FI surface, it introduces a Zeeman field affecting the U 
part of the Dirac Hamiltonian^^— 



2 



rhs (r - t z ) , 



(2) 



where 



= (m x ,m y ,m z ) is FI's magnetization axis, 
s = {s x ,Sy 7 s z ) are the spin Pauli matrices and E z is 
the effective strength of the EPI (absolute magnitude of 
the Zeeman splitting). 

We consider now a central (C) barrier region of length 
Lc, made of BG subject to EPI and described by H = 
Hq + h m + Uo, with Co a possible potential shift, while 
the left side (LS) and the right side (RS) leads are semi- 
infinite normal BG described by Hq. In LS and RS re- 
gions, the dispersion curves are degenerate in the spin 
degree of freedom. In the C region, a spin-splitting E z 
arises between the spin components which are parallel 
{/•) and antiparallel ( k /) to rh, see Fig. [TJ A detailed 
description of the eigenstates of Hq and H is given in 
Appendix [A] 

In the present paper we use E z = 8 meV, close to 
the estimation of Ref. [Tt], and, when not stated other- 
wise, |A| = 150 meV and a temperature of T = 1 K. 
In Fig. [ljb), we show the lowest conduction and valence 
bands in the C region (full, dashed lines) for Uq = and 
in the leads regions (dotted lines). The spin-splitting for 
conduction and valence bands is proportional to the lo- 
calization of their respective states on the U layer (see 
Fig. [5]). Therefore, inverting A will invert the spin split- 
ting for electrons and holes. 

The system is assumed to be invariant under transla- 
tions along Y and the scattering is elastic. Therefore k y 
and E are conserved quantities. We briefly outline here 
the procedure used to calculate the transmission through 
a single barrier by spinor matching and the conductance 
of the system, which has been introduced in Ref. [25J. For 
a given E and k y , solving the Hamiltonian Hq(H), we ob- 
tain an analytical description of the spinors of the prop- 
agating and evanescent modes in the LS and RS (C re- 
gion), which contribute to the scattering state. In the RS 
and LS regions, disregarding the spin, there are four pos- 
sible values of the wavevector k x compatible with a given 
k y and E: k x and — k x , which are propagating modes, 
k x and — k x , which correspond either to propagating or 
to evanescent modes having a finite imaginary part2£. In 
the central part — due to the EPI — the secular equation 
for H leads to spin-dependent solutions of the wavevec- 
tor k x = a n with n = 1, 2 ... 8, described by the spinors 
^ an (x), which are eigenstates of H. 

For an incoming particle of wavevector k x , k y and 
spin-polarization s (a vector describing the up and down 
spin components with respect to the Z axis), we solve 
the linear system determined by imposing the conti- 
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FIG. 2: (Color online) Conductance of spin up and spin down 
carriers as a function of chemical potential fi for a BG device 
acting as a SF calculated for increasing values of Lc- A FI 
is in contact to the U layer of the BG, giving rise to an EPI 
with E z — 8 meV. The gate-bias is A = —150 meV and a 
potential shift of Uo — 5 meV is applied in the C region. In 
the inset, we show the spin-resolved conductance of the SF 
with Lc = 400ao calculated for different temperatures. 



nuity of the scattering state at x = and x = Lc- 
This fixes the output transmission(reflection) coefficients 
tf,ti,tf, t±(rf, r±,r-t,r±) for the allowed k x and k x modes 
in up or down spin orientation. We define the spin- 
resolved transmission probability T\ \> as the sum of 
transmission probabilities in all the outgoing propagat- 
ing modes {\t\\ 2 , \t\\ 2 ), calculated for s compatible with 
a spin-polarization of the incoming particle A'. 

We calculate the conductance of the ballistic system in 
linear response. The 2D two-terminal conductance is 



G\\i — 



(20 



dk x dk y Tx iX 'V x 



df(E - M ) 
dE 



(3) 



with .9 = 2 (accounting for the valley degeneracy), f(E) is 
the Fermi-Dirac distribution, v x the group velocity along 
the transport direction, E the particle energy, and /i the 
electrochemical potential. 



II. SPIN-FILTER 

We now analyze the behavior of the device as a SF. 
In particular, the device acts on unpolarized incoming 
particles, filtering the component antiparallel {^/) to rh. 
In Fig. [21 we show the spin-resolved conductance of the 
device as a function of n, where we choose rh along Z and 
a potential shift of Uq — 5 meV. When ii falls between 
the spin-splitted bands (between 70 and 78 meV) in the 
C region, is exponentially suppressed as a function of 
Lc with an average effective decay length of the order 
of 50 nm, while does not vary. This behavior is due 
to the fact that, in this energy range, transmission of 
spin down particles occurs through evanescent modes, 



9 



05 

CD 
O 
C 

B 3 
§ 10 

"D 



1 


70 


75 






E I 








xo.1 








L„ = 150, 175, 200 (a.) " 
r - 








OFF : 




0.01 










20 25 
H (meV) 


4 y/y — • ***- •' 










J *- *-^ 

N, -» _ *— , 




G tt 










- - G+t 






V 

L c = 200, 175, 150 (a ) - 



E, = 4 meV 



E, = 8 meV 



E, = 16meV 



65 



70 



75 
|i (meV) 



80 



85 



FIG. 3: (Color online) Conductance with transmission into 
spin up Gff and spin down G^ carriers for a BG device acting 
as a SR in ON state, feeded by f-polarized electrons. In the 
inset the spin- flip conductance fraction \ = Git /{Git + G-|--|-) 
is shown for ON (A < 0) and OFF (A > 0) state of the SR for 
Lc = 175ao. The dashed curves, data set (a), are calculated 
with |A| = 40 meV for L c = 200a . 
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FIG. 4: (Color online) In (a), (b) and (c), we present the 
X factor calculated for a SR with different values of the E z 
parameter for Lc varying from 150 to 300 ao- For comparison 
we show in (d), (e) and (f) the corresponding first conduction 
band dispersion: in a dotted curve that of the leads, while in 
full and dashed lines the and ,^-spin-splitted bands in the 
C region of the SR, corresponding to the ON state. 



which exponentially decay (in C). Thus, in the spin- 
splittcd C region and for T -C E z /ks, the device acts 
as an efficient SF(f), i.e. it lets pass only current with 
t spin polarization. Such a SF can be used to generate 
a spin polarized current out of an unpolarized one. Or 
reversely, it can be used as a spin analyzer which detects 
the degree of spin-polarization of charge carriers. This 
possibility will be exploited later in this paper. In the 
inset of Fig. we show the spin-resolved conductance of 
the BG SF for increasing temperature in the range of 1 
to 10 K. As expected, thermal excitations degrade the 
SF efficiency. In particular, the value of E z imposes a 
maximum operating temperature for the device of T ss 
E z /8kB, which corresponds to about 12 K, for E z — 
8 meV. 



III. SPIN-ROTATOR 

We have shown in a previous work2£, that a BG in con- 
tact with a FI can act as an electric-field switchable SR. 
The control of spin- rotation with the gate bias essentially 
depends on the degree of wave function localization on 
one of the BG layers near the Mexican hat energy dis- 
persion region. A useful parameter to characterize spin- 
rotation is the ratio x of the conductance associated with 
a spin-flipped transmission to the total conductance^. 
As shown in the inset of Fig. [3J we can put the device 
OFF (xoff « 6%) or turn it ON (xon ~ 60%) by revers- 
ing the gate bias. The performance of the spin-rotator 
is basically limited by the finite \ fraction in the OFF 
state, due to the non perfect layer localization of elec- 
trons contributing to transport. With a smaller gate-bias 
of |A| = 40 meV— , the layer confinement is more effec- 



tive (see Fig. [5]) and our model predicts a moderate per- 
formance enhancement (xoff ~ 4% and xon ~ 80%). 
Fig. [3J shows the spin resolved conductance in the ON 
state for a SR feeded by a f-spin-polarized source lead, 
for increasing lengths of the magnetic barrier. In par- 
ticular, for L — 175ao (ao = 1,A2A the in-plane nearest 
neighbor distance) , a strong spin-flip resonance is present 
at /i w 74 meV, which is well inside the SF operational 
regime presented in Fig. [5] 

The origin of the spin rotation is easily explained. 
A f-spin-polarized electron is described inside the spin- 
rotator as the superposition of two components with spin 
polarization and */ (polarizations which are parallel 
and antiparallel to the FI magnetic axes and are eigen- 
states of the EPI system). These two components, being 
coupled differently to the EPI, travel with a different k x 
wavevector and accumulate a phase difference Ak x ■ Lc 
in a single crossing of the C region. The phase differ- 
ence translates into a net rotation of the initially f-spin- 
polarized electron. Due to the complex 2D BG dispersion 
curve, it is difficult to establish an immediate relation for 
the spin-flip resonance condition between the parameters 
A, E z and Lc- However, for a wide range of E z values, 
a spin-flip resonance is observed in the Mexican hat re- 
gion of the ^-spin-polarized band, for |A| = 150 meV 
and L c « 200ao, as shown in Fig. @] The order of mag- 
nitude of Lc is related to 7r/fc m i„, where k m i n is the 
wavevector corresponding to the minimum of the first 
conduction band (Eq. I A10|) . In Fig. [H we compare the 
Xon and xoff of three systems with Ez = 4 (a), 8 
(b) and 16 meV (c), for a C region of length Lq = 150, 
200, 250 and 300 o . We also plot in Fig. gjd-f) the 
corresponding k y — dispersion curves for the normal 
BG (dotted) and for the BG subject to EPI in the ON 
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state, in a full curve for /~- and in a dashed curve for 
i/-spin-polarization. In correspondence to the edge of 
the ^-spin-polarized band, transport is often associated 
with a maximum of xon, corresponding to the fact that 
the majority of the electrons which tunnel through the C 
region are spin- flipped (i.e. they satisfy Ak x ■ Lc ~ 7r). 



IV. SPIN FET 

We now calculate the total conductance G = 
Ea A' Caa' for a hybrid setup made by the series of a 
SF(t), a SR(ON/OFF) and a SF(|), each one built from 
BG in contact with a FI with fh along Z, Y and — Z, re- 
spectively. Ideally, the SF(f) selects the f-component of 
the incoming unpolarized electrons, resulting in a spin- 
polarized current. The SR introduces a spin precession 
which we can turn ON or OFF with the gate bias (see 
inset of Fig. [3] and related discussion). Finally, the SF(4_) 
measures the degree of spin-rotation, because it (ideally) 
lets pass only carriers which have been spin-flipped by 
the SR. Therefore this structure realizes a complete spin- 
tronic scheme of creation, manipulation and measure- 
ment of spin-polarized currents, which does not require 
spin-polarized leads. The calculation of the transmission 
is performed by applying the transfer-matrix (TM) for- 
malism, which we briefly review in Appendix[Bj For each 
barrier we obtain the corresponding TM, which requires f 
and t for the individual scattering problems of a particle, 
approaching the barrier from the LS or from the RS. The 
TM links the LS modes to the RS modes and therefore 
is multiplicative, in the sense that the TM of a series of 
barriers is the ordered product of the corresponding indi- 
vidual TMs. We thereby neglect the contributions of the 
evanescent modes connecting different scattering regions 
of the system. Therefore we consider the three spintron- 
ics blocks separated by a distance of L = lOOOao, where 
any effects of such evanescent modes are negligible 3 - 7 .. 

Fig. E] shows the total conductance of a spintronic de- 
vice made by the composition of SFs and a SR as a func- 
tion of ijl. The full black lines, as indicated, show the con- 
ductances of an "open" series, i.e. the SF(f)-SF(t), and 
of a "closed" one SF(t)-SF(i). For fj, in the operational 
region of the SF, the total conductance of the "closed" 
series is suppressed by almost 4 orders of magnitude with 
respect to the "open" one. The remaining lines represent 
the conductance of a "closed" series including the SR de- 
scribed in Fig. [3l The dashed green line — the SR(ON) 
case — is quite close to the conductance of the "open" se- 
ries. The dotted red line — the SR(OFF) case — exhibits 
a conductance suppression by approximately a factor 10, 
with respect to the ON case. This expresses a measure 
for the efficiency of the SR, for which, in fact, we note 
that Xon/xoff ~ 10 for the corresponding data set (b) 
in the inset of Fig. [3J 
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FIG. 5: (Color online) Hybrid devices combining the SF and 
SR blocks described in Figs. [2] and O respectively. The full 
lines represent the total conductance for a series of two SFs 
in "open" or "closed" configuration. The dashed and dotted 
lines are for the total conductance of combined system made 
by the series of a SF(t), a SR(ON/OFF) and a SF(|). The cal- 
culation has been performed for a gate-bias of | A| = 150 meV, 
with L c = 400a for the SFs and L c = 175<z for the SR. 



V. DISCUSSION AND CONCLUSIONS 

We now discuss the possibility of actually realizing the 
EPI in graphene devices. One of the few concrete ex- 
amples of a FI is EuO. The first realization of the EPI 
coupling, originating a spin-splitting, has been experi- 
mentally proven in a EuO/supcrconductor interface 27 . 
The possibility to incorporate FI in nanostructures has 
been recognized to be extremely attractive for the real- 
ization of spintronic nanodevices and recently much ef- 
fort has been put into the development of the FI tech- 
nology. In particular, important steps have been made 
in the control of the epitaxial growth of EuO on Si and 
GaAs— £SL. EuO seems promising with its semiconduct- 
ing gap of about 0.7 e V 30 i 31 and the possibility to be 
grown in thin films of a few nm thicknes o 31 ' 32 . Regard- 
ing the practical realization of the device, a suitable FI 
should have a sufficiently large bandgap and retain its 
properties when grown in thin films. 

The occurrence of a EPI for graphene deposited on 
a FI has been proposed by several authors^—, with a 
tentative estimation of the expected Zeeman splitting of 
E z s» 5 meV— . When a BG is placed in contact with 
a FI, the EPI, being short-ranged, affects only the U 
graphene layer, which is in direct contact with the Flia. 
In fact, applying the contact exchange model between 
magnetic ions and itinerant electrons proposed in Ref. HH, 
and using the asymptotic atomic wave functions for car- 
bon 3 ^, it is easy to show (see Appendix [Cj) that the ratio 
between the EPI strength on the L and U layer is on the 
order of e~ 2Kd ° ~ 10~ 3 , with d — 3.4 A the interlayer 
distance in BG and k sa O^la^ 1 the asymptotic exponent 
lorG 3 !. 
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In summary, we have demonstrated that the exchange 
proximity interaction in bilayer graphene in contact to 
a ferromagnetic insulator can be exploited as a means 
for electrical spin manipulation. We have shown that 
this system acts both as a switchable spin-filter or spin- 
rotator, which are basic building blocks for spintronics. 
As an example, we have shown how to realize a complete 
spintronic structure for the creation, manipulation and 
detection of spin currents — a spin FET, out of an ini- 
tially unpolarized stream of electrons and calculated its 
operational efficiency with a transfer matrix approach. 
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1. Spinors of BG without EPI 

In this section, we give the analytical expressions for 
the spinors of normal BG, i.e. without EPI and neglect- 
ing the trigonal warping effects. Spin is degenerate and 
therefore we can consider xx as scalar complex numbers. 
The Hamiltonian system (Ho — E)*f> = leads to the sec- 
ular equation 



t' 2 (e' 2 



5 Z 



0(A5) 



with k 2 = k 2 + k 2 . If we solve for the energy we obtain 
the BG eigenstates 



(4) 5 



t' 2 



X2U2 

l/ + Y ±t 'i- + 4 — + k2 - (A6) 



If, instead, we solve for k x we obtain the BG modes con- 
sistent with energy E and transverse wavevector k y 



Appendix A: Bilayer graphene eigenspinors 



We consider the BG Hamiltonian 



' 2 



H 



V 



Vfk- 

t± 




v f k+ 
4 + h Tl 





t± 


A 
2 



+ 






Vfk- 
A 
2 



(Al) 



with k± = k x ±ik y and where A, vt and t± have been in- 
troduced in SectionQ] EPI affects only the U plane and is 
contained in h m (Eq. ^ , all other terms are proportional 
to the identity in the spin subspace. The Hamiltonian 
acts on the spinor 



* = 



J/q I ^ik y y 



Xb> 

XA' 



\J L X L 



(A2) 



y 



where A, B refer to the two inequivalent sublattices on 
the U BG layer, A' ,B' to that of the L one. L x and L y are 
the channel dimensions along X and Y directions. Now 
we distinguish the two spin components along the Z axis, 
perpendicular to the plane, therefore \x, with X = A, 
B, A' , B' , has to be regarded as a two-component spinor 



-k 2 + e' 2 + 6 2 ±t\ e /2 -S 2 +4 



e' 2 S 2 



(A7) 



The spinor components can be expressed (for e' ±5) 



as 



XA 
XB 
XA' 
XB' 



XB' 



t'(e' + 5) 

t'{e' 2 - S 2 
(k x + iky) 
e' + S 
t'{e' 2 



-XB' 



XB' 



y/A\B\ 2 + i' 2 (e' -5) 2 C' 

where we adopted the following notation 

A = (e 1 -5) 2 + \k x -iky\ 2 

B = (e' + S) 2 -k 2 

C = [e 1 + 8) 2 + \k x +ik y \ 2 



V 



(e' - S) 2 - fc 2 . 



(A8) 



a. Layer localization and trigonal warping corrections 



Xx 



( 4>xt 



(A3) 



We introduce the following notation, similar to that 
chosen in Ref. [2a, 



The probability to find the electron, of a specific eigen- 
function, on the U plane is given by 



Ixa| 2 + Ixb| 2 



i 



1 4- .. -4B 2 .. 
1 + t'^(e'-8)^C 



(A9) 



U — r>i 



Vv-Vl 



2hvr. 



E 

hup 



"° — 2hv F 



(A4) 



and correspondingly P L = 1 — P u . For our purpose, 
the most important part of the dispersion curve is the 
Mexican hat region (|fc| w fc m m). For this reason we 
now describe in more details this minimum for the first 
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FIG. 6: Spinor-projection on the U BG plane for the eigen- 
state corresponding to k m i n of the first conduction (full line) 
and valence band (dashed line), as a function of the layer 
potential bias A. 



conduction band. Its wavevector, energy and layer pro- 
jection are given by 



P. 



rain 
U 



±SX 



X 2 



i , t' 2 +4S 2 1+X 2 +X ( X 2 iV 
1H IP— 1+X 2 -X l A ~ 1 ) 



(A10) 



with X = ^ iS \ +t ,i ■ In Fig. [6j we show the layer local- 
ization properties of the first conduction band minimum 
with a full curve (and the first valence band maximum 
with a dashed line). For a sufficiently small |A|, the cor- 
responding states are strongly localized on the upper or 
on the lower graphene plane, depending on the sign of 
the applied bias A. Similarly, the low energy states of 
the first conduction band share analogously strong layer 
localization properties with the conduction band mini- 
mum. In particular, the layer localization complete for 
the k = state (P u = 1 or 0). 

We analyze, now, the effect of the trigonal warping 
correction^ on the layer localization properties of BG. 
The trigonal warping correction, acting on the spinor in 
Eq. EH is 



H 3 = v 3 



Co 



■ {k x 



kyTy), 



(AH) 



with er, r the sublattice and layer Pauli matrices, respec- 
tively. We have solved the system with a finite trigonal 
warping with v 3 /vf = 0.1—, found its eigenstates and 
calculated their layer projections. In Fig. [3 we compare 
the projection on the U plane for the system without 
trigonal warping (solid curve) and including the trigonal 
correction (dashed curve) for the dispersion along the X 
axis in (a) and Y axes in (b) for A = 40 meV, energy 
for which the isotropic Mexican hat dispersion is heav- 
ily distorted to a trigonal symmetry. In both cases, the 
low energy states are essentially localized on the L layer. 
We conclude that the layer projection properties are only 
slightly affected by the presence of the trigonal warping 
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FIG. 7: Projection of the eigenstates of the first conduction 
band on the U BG plane, with (dashed line) and without 
(full line) trigonal warping corrections, for an applied bias of 
A = 40 meV. (a) is along the X axis and (b) along the Y 
axis. 



in the BG Hamiltonian and the mechanism which permit 
to turn ON or OFF the BG spintronic functionalities is 
still valid. 

The presence of the trigonal warping can however ef- 
fect the transmission of the system. The spin filtering 
properties are unaffected, as long as a gap is present 
with spin-splitted bands. The spin rotation effect is due 
to the interference effect between different spin compo- 
nents propagating with different wavevectors from the 
Mexican hat spectral region. The trigonal warping dis- 
tortion of the Mexican hat can change the accumulated 
phase-difference and therefore leads to differences in the 
form of the spin-flip transmission resonance. The dis- 
tortion is more pronounced at smaller applied gate bias 
| A [22, inducing a progressive angle dependence of the 
spin-rotating properties. In principle, also in this small- 
gap regime, spin rotation could still be exploited in 
angle-selective transport experiments. Instead, for suf- 
ficiently large gate bias (and therefore semiconducting 
gap), the trigonal warping effect is limited to a minor 
distortion of the Mexican hat and the properties of the 
system are essentially unchanged. This scenario is met 
for |A| = 150 meV used in the paper. 



2. Spinors of BG with EPI 

In this section we derive the analytical expressions for 
the BG spinors in the case of a finite EPI acting on the 



7 



U plane. EPI is not diagonal in the spin components and 
the spin variables are individually addressed. Still the 
Hamiltonian system (H — E)"^ — 0, can be analytically 
solved leading to the secular equation 



and RS regions, for modes which are incoming towards 
the scattering center and outgoing from it, respectively. 

We define the scattering matrix S and the transfer ma- 
trix T through the following relations 



'2\-| 2 



= — [2B( S '-S)-t' 2 (e' + 6)] 



The secular equation solved for k x has in general 8 com- 
plex solutions 



-Pi 



with 



+ 
\ + J 



Pi 



P2 



\ + J 



PI 



P 3 



f-_\ 

+ 
\ + J 

(A12) 



B-V- — 
4 



P 2 = ^{e' 2 -5 2 ) 



P 3 
Pi 



-BV-t l2 {e' 2 



a 



-B 



a 
~2 



2B{e' -8)-t' 2 {e' + 5) 



where B = B + k 2 and V = V + k 2 . Real solutions 
correspond to propagating modes in the region with EPI 
interaction while modes with an finite imaginary part 
give exponentially decaying modes at the border of the 
EPI region. 

We omit the the expressions for the spinor components. 
We obtained them with straightforward derivation, from 
the secular equation, as a function of <pB'±, which in a 
second time we numerically fixed using the normalization 
condition. 



Appendix B: Multiple barriers: 
method 



transfer matrix 



Let us consider a ID channel with a finite number of 
modes M in which a scattering region is present. For each 
scattering center, we can fictitiously divide the system in 
a left-side (LS) and a right-side (RS) leads, which we 
assume to be semi-infinite. The wavefunctions in the LS 
and RS leads are described by 



(Bl) 
(B2) 



a 



where 



is the o-th mode of the channel carrying an 
unity of current, where ± stands for forward-going (+), 
i.e. from the LS to the RS, and backward-going (-). 
j{l,r) anc j q(l,r) are £ ne coe fg c i en t vectors, in the LS 



s ( liL) 

f(L) 



(B3) 



It is immediate to identify the elements of the scattering 
matrix with reflection and transmission coefficients, so 
that 



t r' 
r t' 



(B4) 



but we are interested in the transfer matrix because it 
is multiplicative, i.e. the transfer matrix of a series of 
scatterers is given by the ordered multiplication of the 
individual transfer matrices, for each one of the scattering 
centers, or 



(B5) 



We can obtain an expression for the transfer matrix in 
terms of reflection and transmission coefficients by com- 
paring the action of the S and T matrices in Eq. 



q(L) = t ,f(R.) + r f(L) 

m = To 3 dW+T 2 ,ijW 



(B6) 



We obtain 



T2,2 

T 2 ,i 

T l,2 



[fr 1 

t-r'[t'] _1 r 
r'[tr X - 



(B7) 



In practice, in order to calculate the transfer matrix for 
each one of the scatterers, the transmission and reflection 
coefficients for a particle approaching from the LS and 
RS of the scattering center are needed. 

From the total transfer matrix, it is then possible to 
obtain the transmission and the reflection properties of 
the overall system by the following relations 

t' = [T2.2]" 1 

r = -[T 2 ,2] _1 T 2 ,i 
t = Tij — Ti : 2 [T2,2] 1 T24 
r' = T ll2 [T 2 ,2] _1 . (B8) 
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1. Properties of S and T matrices 

The scattering matrix has to be unitary in order to 
ensure charge conservation in a barrier, i.e. Jj n = I ou t- 
Explicitly 



( T T V st s ( fiR) 



which is satisfied if S^S = 1 and therefore = [S] _1 . A 
well-known consequence of the unitarity of the scattering 
matrix is that |t| = |t'| and |r| = |r'|, valid for any kind 
of elastic scatterer. 

The condition for a stationary equilibrium current 
through the barrier is given by 



valence bands can essentially be described by the carbon 
atoms p z -orbital. The tails of the graphene p z orbital en- 
ter the FI, where the magnetic ions are distributed. We 
assume an homogeneous distribution of magnetic ions in- 
side the FI, with density rii ons in the FI with fixed average 
spin polarization (Si) along Z. We obtain the following 
exchange potential for the graphene electrons 



V ex {r) = -jn ions Q(z - L)S z (Si), 



(C2) 



where <d(x) is the step function and z = L identify the 
surface of the FI. 

An established result in atomic physics 3 ^ is that we 
can describe the asymptotic behavior of wave functions 
for valence electrons in an atom, at large distances, as 
tp(r) = Rn t m(r)Y m! i, with R(r) 



I = 



jXR) - _ Q{R) 



(B9) 



R(r) = Ar 1 '^ 



(C3) 



which imposes the following property on the transfer ma- 
trix T^cr z T = a z . 



with k = V 21, where / is the ionization potential for that 
electron in the atom. A p z -electron on a graphene plane 
will therefore be described as 



Appendix C: Estimate of EPI effects on the distant 
BG layer 

The goal of this section is to provide an estimation of 
the relative importance of the EPI of the two layers of 
BG, when the U layer is placed in direct contact to the 
FI surface. We will consider the graphene layers and the 
FI surface oriented along the XY plane with the FI sur- 
face at z = 0, the U plane centered at z — Ljj and the 
L plane centered at z = Ll, with Ll — Ljj — dg the BG 
interlayer distance. As noted in Ref. [HI, the exchange 
coupling between an itinerant electron and the local mo- 
ments in FIs (like EuO) typically dominates the coupling 
to the magnetization. The exchange potential for a mo- 
bile electron, arising because of the exchange interaction 
with core electrons in a magnetic ion, is modeled by 33 



ions 



V ex (r) = - ] ^5(\r-R i \)S i S, 



(CI) 



where Si is the ion's total spin, S is the electron spin op- 
erator and j an exchange energy parameter. This expres- 
sion is used in Ref. [33| to model the exchange potential of 
conduction electrons due to the presence of localized core 
d-electrons in Mn magnetic ions. The assumption is that 
the wavefunctions of mobile electrons can be considered 
approximatively constant in the range of variation of the 
magnetic ion's occupied orbitals. 

In our model the mobile electrons are the graphene 
bilayer conduction and valence band electrons. In the di- 
rection perpendicular to the graphene plane, this mate- 
rial is practically one-atom thick and its conduction and 



ip(r, 6) oc r 



l/«-l e -r« cos( 



(C4) 



The direct EPI energy between the FI and a graphene 
7r-orbital from the U graphene plane is proportional to 



E ( P oc 



■ 11k -2rn 

dr r ' e 

L-L-i -In 



d6 cos 



Xi, 



with t — r/(L — Ljj) and 



Xi = 



2 - 1 1 

-r h arccos — 

t t 



(C5) 



(C6) 



? { P, where Ef ] is the 



We now calculate the ratio E z /E 
magnitude of the EPI with a carbon p z orbital from the L 
graphene plane, which is further away from the FI surface, 
than the U one. Observing the form of the EPI in Eq. IC5I 
and using the fact that L — Ll > L — Lu, we can infer the 
following condition for E 



(L) 



E 



and therefore the ratio 



-2(L U -L t ), 



E 



(07) 



(08) 



For two neighboring graphene layers (interlayer distance do 
around 0.34 nm), and employing the value for C atoms ft = 
0.910 a^ 1 from Ref. H 

we conclude that the ratio of the 
exchange interaction (Jl/Ju) is of the order of 2 x 10 -3 . We 
can therefore safely neglect the EPI effect on the lower layer. 
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